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Abstract 

Why reaction rate constants for enzymatic reactions are typically inversely propor- 
tional to fractional power exponents of solvent viscosity remains to be already a 
thirty years old puzzle. Available interpretations of the phenomenon have not led 
to consensus among researches about its origin. They invoke to either a modifica- 
tion of 1. the conventional Kramers' theory or that of 2. the Stokes law. We show 
that there is an alternative interpretation of the phenomenon at which neither of 
these modifications is in fact indispensable. Basing on an analogy from the theory 
of adsorption on heterogeneous surfaces we reconcile 1. and 2. with the experimen- 
tally observable dependence. We assume that an enzyme solution in solvent with 
or without cosolvent molecules is an ensemble of samples with different values of 
the viscosity for the movement of the system along the reaction coordinate. We as- 
sume that this viscosity consists of the contribution with the weight q from cosolvent 
molecules and that with the weight 1 — q from protein matrix and solvent molecules. 
We introduce heterogeneity in our system with the help of a distribution over the 
weight q. The function of the distribution is a unique characteristic of the solution 
of enzymes (type of the enzyme, cosolvent molecular weight, pH, temperature, etc.). 
We verify the obtained solution of the integral equation for the unknown function of 
the distribution by direct substitution. We conclude that even at linear relationship 
between the solvent viscosity and that for the movement of the system along the 
reaction coordinate our approach enables us to obtain the required dependence. All 
parameters of the model are related to experimentally observable values. The mean- 
ing of fractional exponents appears to be the characteristic for the behavior of the 
distribution with the variation of the weight q. Our approach yields the existence 
of the limit value for the fractional power exponent with the decrease of cosolvent 
molecular weight that is in agreement experimental data known from the literature. 
This limit value is determined by the properties of the protein structure and thus 
is a unique characteristic of the type of the enzyme only rather than that of the 
solution. General formalism is exemplified by the analysis of literature experimental 
data for oxygen escape from hemerythin. 

Key words: enzyme catalysis, solvent viscosity, Kramers' theory, protein 
dynamics. 
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1 Introduction 



Viscosity dependence of enzymatic and protein (ligand binding/rebinding) 
reactions is a long standing problem for biophysics pQ, [2], [3J, [I], [5], [6], [7j, 
0, H, [10], [EE], [12], [13], [H], [15], [H>], [IT], PS], [H]. For such reactions 
the functional dependence of the reaction rate constant for the rate limiting 
stage k on solvent viscosity rj has the form 

k OC g (1) 

(v/Vo) 



where t]q is the viscosity of pure solvent (for water i] = 1 cP at room tempera- 
ture) and < (3 < 1 (usually (3 ~ 0.4-h0.8). This dependence is experimentally 
verified in the range of variation of solvent viscosity by two orders of magni- 
tude ?7 < 100 cP. Similar dependence also takes place for folding of proteins 
(see [20], [21], [22] and refs. therein) and at formation of protein structure [23] . 
However we will not touch upon these processes in the present paper. 

Much efforts were devoted to explaining the functional dependence ([1]) in the 
previous century with no obvious consensus and commonly accepted mecha- 
nism of the phenomenon. For instance the authors of [13] conclude that "there 
seems to be no general agreement yet about the origin of the fractional (3 value 
in Eq.l" . The authors of [TB] draw to a similar conclusion: "At present there 
is no general agreement on the meaning of fractional exponents implying a 
breakdown of Stokes law" . Despite steady growth of experimental data little 
has changed in this issue since the date of the cited papers (cf., e.g.,: "The 
meaning of fractional exponents implying that the friction coefficient does not 
vary linearly with the solvent viscosity, essentially a violation of Stokes law 
(/ oc r] s ), is not clear at present" [22]). Detailed studies revealed that in fact 
the fractional index of a power (3 is a function of cosolvent molecular weight M 
(i.e., the mass of a cosolvent molecule expessed in atomic units and measured 
in Daltons) (3 = (3(M) [13]. If one varies the solvent viscosity by large cosolvent 
molecules with high molecular weight that do not penetrate into enzyme then 
one obtains that the fractional exponent (3 — > 0, i.e., the reaction rate constant 
does not depend on solvent viscosity. With the decrease of cosolvent molec- 
ular weight the fractional exponent (3 increases. In the limit of hypothetical 
"ideal" cosolvent with very small molecular weight (cosolvent molecules freely 



Email address: sitnitsky@mail.knc.ru (A.E. Sitnitsky). 



2 



penetrate into enzyme and are distributed there homogeneously) it tends to 
some limit value 

(3 max = lira (3{M) (2) 

M -tM m ;„ 



For oxygen escape from hemerythin this value is (3 max ~ 0.79. The latter is 
neither experimental value nor a calculated one. It is an extrapolated number 
(see [13J for details). The existence of enzymatic reactions with (3 = 1 |17j . 
[2T] suggests that the limit value is a unique characteristic of the type of the 
enzyme and can take the values up to 1 as well. 

The main tool to describe the viscosity dependence of a reaction rate constant 
is the high friction limit (also called strong damping or overdamped regime) of 
the Kramers' theory [21]. In it the reaction is conceived as a diffusion process 
of a particle with some effective mass along a reaction coordinate over some 
potential surface. The friction coefficient for the particle v is supposed to obey 
the Stokes law v oc r\ that yields the well known dependence 

k K oc —j— (3) 
V/Vo 



As was stated above for enzymatic reaction this dependence is inconsistent 
with observations. There are two strategies to interpret the experimentally 
observed dependence ([T]) at present: 1. to modify the Kramers' theory and 
2. to modify the Stokes law. The first approach usually leads to rather com- 
plicated theoretical constructions that have not proved to yield a universal 
and commonly accepted resolution of the problem. It is realized in, e.g., the 
Grote-Hynes theory (GH) [23], the Zwanzig model [26J or the model suggested 
in [27], [12]. GH gives that the rate dependence on solvent viscosity should 
be weaker than that predicted by Kramers' one. The authors of this model 
argue that the friction coefficient should be proportional to a high-frequency 
viscosity at an appropriately renormalized frequency [25]. Assuming a pos- 
itive power dependence of viscosity on frequency, they find that the renor- 
malized high frequency viscosity has a fractional power dependence on the 
low-frequency viscosity which can be measured as the usual viscosity rj. Thus 
their assumption is in fact introducing the fractional exponent "by hand" in 
a phenomenological manner. The Zwanzig model infers the reaction rate con- 
stant from the first principles but yields too small value for the fractional 
exponent (3 = 0.5. The model [27], p2| deals only with the limiting case of the 
hypothetical "ideal" cosolvent with very small molecular weight and yields for 
the fractional exponent the limit value (3 max ~ 0.75 that is in good agreement 
with the extrapolated one f3 max ~ 0.79 from the paper [13]. No extension of 
this model to the case of realistic cosolvent with finite molecular weight has 
been suggested yet. The second approach takes into account that the value of 
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viscosity in enzyme active site differs from that of the solvent [2] and may be 
a nonlinear function of the latter. It is suggested that the fractional exponent 
/3 is the degree with which solvent viscosity is coupled with (frequency depen- 
dent friction) [5], [2Hj, [IE], or penetrates into (position dependent friction) 
[I], [E] the protein interior. Any of these variants yields the modification of 
the Stokes law v oc rf . However the fractional exponent (3 in these approaches 
still appears as an empirical parameter. It accounts for the nonlinearity of the 
relationship between the values of the solvent viscosity and that in the enzyme 
active site but its value is not calculated and its origin remains mysterious. 

In our opinion there is another way to interpret the experimental data that re- 
quires neither modification of the Kramers' model nor that of the Stokes law. 
Thus the aim of the present paper is to reconcile the conventional Kramers' 
model and the Stokes law with the experimentally observed dependence (PQ). 
The main premise of our approach is in the fact that a realistic enzyme solu- 
tion in solvent with or without cosolvent molecules is actually an ensemble of 
enzymes with different conditions (values of viscosity for the movement of the 
system along the reaction coordinate). We show that even at linear relation- 
ship between solvent viscosity and that for the movement of the system along 
the reaction coordinate we can obtain (TjQ) if we take into account heterogeneity 
of conditions in the ensemble. It should be stressed that the present approach 
can be generalized to be based not on the Kramers' model but on the mod- 
els [25] . [26\ or |27j . |19j . However for the sake of definiteness we choose the 
Kramers' model and reconcile the simple dependence (j3J with the experimen- 
tally observable expression for the reaction rate constant (pQ). The latter is 
obtained by averaging of individual Kramers' rate constants over the distri- 
bution. The aim of the present paper is to show that the idea of heterogeneity 
enables one to resolve the thirty years old puzzle of solvent viscosity effect 
on enzymatic reaction rate constant in a conceptually much more simple way 
than modification of either the Kramers' theory or that of the Stokes law. 
Also we stress that our distribution has nothing to do with the nonthermally 
equilibrium form of the reactant distribution during the reaction employed in 
the approach of [29] . 

There is a noteworthy analogy between the problem under consideration and 
that of adsorption on heterogeneous surfaces. In the latter one reconciles the 
Langmuir isotherm grounded by statistical mechanics and giving the num- 
ber of adsorbed molecules 9 proportional to pressure P in the low pressure 
limit 9 oc P with the experimentally observable phenomenological Freundlich 
isotherm giving the dependence 9 oc P 5 where < 5 < 1. One attains this aim 
via introducing a distribution over energy due to heterogeneity of the surface. 
As is well known the surface heterogeneity in adsorption theory is taken into 
account by integral equation approach (see, e.g., [30], [3J] and refs. therein). 
This line in adsorption theory has a long history. It was initiated by the paper 
of Zel'dovich [32] (see, e.g., [33J for authoritative discussion of priority ques- 
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tions). The next important steps were made by Sips [31], [35] and Misra [36] 
who used Stieltjes transform. Detailed discussion of related problems is given 
in [37], [3S], |39j . For the solution of the corresponding integral equations the 
authors of these papers used the so called condensation approximation [37], 
[3T] . Finally this line was culminated by the paper of Landman and Montroll 
[4Ti] who applied powerful Wiener- Hopf technique for the solution of the corre- 
sponding integral equations. The enormous preceding experience accumulated 
in this field of chemical physics provides an invaluable source of information 
for our problem both from the side of conceptual aspects and that of mathe- 
matical technique. 

The paper is organized as follows. In Sec. 2 the character of the relation- 
ship between the viscosity for the movement of the system along the reaction 
coordinate and that of solvent viscosity is discussed. In Sec. 3 we formulate 
the integral equation for taking into account the heterogeneity of conditions 
in the enzymes from the ensemble via introducing the distribution over the 
weight q. In Sec. 4 we cast this equation to the form solvable with the help of 
Fourier transform and verify our solution by direct substitution. In Sec. 5 we 
obtain the relationship of the all parameters of the model with experimentally 
observable values. In Sec. 6 the results are discussed and the conclusions are 
summarized. In Appendix some known mathematical formulas are collected 
for convenience. 



2 Viscosity for the movement of the system along the reaction 
coordinate 

The reaction coordinate in the Kramers' theory is a notion that enables one to 
reduce very complicated dynamics of the system in a multidimensional config- 
urational space to the movement of a particle with some effective mass along 
one effective dimension. For the case of an enzymatic reaction the movement 
of the system along the reaction coordinate is determined by: processes in the 
substrate molecule (e.g., stretching of the bond to be cleaved); friction for 
these processes by the environment in the active site of the enzyme (that is 
stipulated by solvent and cosolvent molecules and also protein side chains lo- 
cated there); relevant motion for catalysis of side chains and larger fragments 
of structure from protein interior (such as, e.g., a-helicies and /3-sheets) and 
friction to it by solvent and cosolvent molecules and also by protein matrix 
in their environment; etc. We denote the viscosity for the movement of the 
system along the reaction coordinate as ( and suppose the conventional Stokes 
law for the friction coefficient of the particle v oc (. The solvent viscosity we 
denote as earlier 77. The linear relationship between 77 and ( can be motivated 
as follows. 
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There seems to be consensus among researches that the effect of solvent vis- 
cosity on the enzymatic reaction is mediated by protein dynamics [2], [3], [5], 
IS] , [H] , [H] , PH] , [EE] , [12] , [15] • That is why we further identify the viscosity for 
the movement of the system along the reaction coordinate with the viscosity 
at relevant motion for catalysis of individual fragments in protein interior and 
in enzyme active site. Experiment yields that protein dynamics over a wide 
range of temperatures obey the relationship between solvent viscosity 77 and 
the mean square displacement (MSD) < u 2 > [41], [42] 



In (77/770) 



where g is a constant. We notice that firstly (JH) is divergent at 77/770 = 1 
and secondly nobody verifies ([1]) experimentally in the limit 77/770 — * 1+. In 
practice one starts from some values noticeable larger than 1. For instance in 
[T3] {r}/Vo)min ~ 1-6. That is why further on we consider the range of solvent 
viscosity 



The theoretical counterpart of (j4J) given by, e.g., the most simple and extraor- 
dinary successful "model of bounded diffusion" [43], [44], [5] yields for the 
MSD the expression 

< u 2 >= 6 
mujn 



where is the Boltzman constant, T is the temperature, m is the mass of 
the moving fragment and u is the frequency defined by the curvature of the 
potential well near its bottom for the motion of the fragment. This motion 
is determined by the value of the local viscosity ( of the environment for the 
motion of the fragment in protein interior. 

First we notice that a nonlinear relationship between ( and 77 is compatible 
with (j4J) and ([6]) with the linear one being a particular case. We assume that 
C/i]o consists of the contribution with the weight 1 — q determined by protein 
matrix and solvent molecules and the relative contribution with the weight q 
from cosolvent molecules in the fragment environment. The weight q may take 
values from the range 

< q < 1 (7) 



The term "relative" means that the latter contribution must tend to zero 
for the case of pure solvent (without cosolvent molecules). That is why for 
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it we adopt the conventional Andrade-Eyring-Frenkel' form exp \Ej (/c#T)] 
subtracted by 1 (because for the case of pure solvent 77/7/0 = 1 or equivalently 
at E = the exponent is 1). Here E is the activation energy (E > 0). Thus 
we consider the form 



C/Vo = a\(l-q) + q 



exp 



E 



(8) 



The parameter a is related to the fractional power exponent (3 that will be 
obtained later on theoretical grounds. Also below the explicit dependence 
a = a (M) of the parameter on the cosolvent molecular weight M will be 
obtained from the analysis of experimental data. 

One can easily see that the nonlinear relationship between the solvent viscosity 
7] and the local viscosity for the motion of the fragment in protein interior ( 
of the form 



C. 




q)+q 




- 1 



(9) 



is compatible with the relationships (jlj), (jSj) and © at any k > if we set 
E 



muo^K 



(10) 



At K = 1 we have the particular linear case 

- = a[(i-g) + gf--i) 



There in no wonder that we can obtain the experimentally observable depen- 
dence ([I]) for the nonlinear case < k < 1 simply by setting (3 = K. Our aim is 
to stress that even the linear case k — 1 can lead to the required dependence 
([1]) as well. At the same time it should be emphasized that the assumption 
about linearity is by no means crucial for the validity of our approach. The 
latter can still remain workable for the general nonlinear case even at k 7^ (3. 
However we will not explore this option in the present paper and restrict our- 
selves by the linear 1. In our opinion namely this case can elucidate 
the essence of our approach with full clarity. We pursue this line in the next 
Sec. 
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3 Formulation of the integral equation 



The value of viscosity for the movement of the system along the reaction 
coordinate may be different for different enzyme molecules (samples) from the 
ensemble of enzymes. In the individual samples from the ensemble we have 
heterogeneity of conditions both in the enzyme active sites and in protein 
interior. That is why the experimentally observable dependence of the reaction 
rate constant on solvent viscosity should be determined by some averaged value 
of the viscosity for the movement of the system along the reaction coordinate. 
We consider the following way to take into account this heterogeneity by 
averaging with the help of some integral equation. 

We assume that due to heterogeneity of conditions in the individual samples 
from the ensemble of enzymes the weight q in (jSj) and consequently in ffTTT) 
can take different values from the range < q < 1. Thus the ensemble is 
characterized by the distribution p(q) over the values of the weight q. We also 
assume that the reaction rate constant for the individual sample from the en- 
semble of enzymes firstly obeys the Stokes law v oc C/^o an d secondly is given 
by the Kramers' formula Q that is oc (C/Vo)~ ■ Then the experimentally 
observable reaction rate constant is the average over the distribution 




The distribution p(q) must be normalized 
i 

/ dq p(q) = 1 (13) 



As we want the reaction rate constant to be the experimentally observable (pQ) 
and the local viscosity ( to be fill I) then the relationship (1121) takes the form 

l Uq- ^ - - (14) 



{V/Vof a { (l-q) + q(v/Vo-l) 



The latter is a Fredholm integral equation of the first kind for the unknown 
function of the distribution p(q). It is analyzed in the next Sec. 
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4 Analysis of the integral equation ( 1141) 



We introduce the critical value of the weight q 

Qc = \ (15) 



Equation ffT4l) can be cast into the form solvable with the help of Fourier 
transform. Making use of N2.2.6.15. from [15] and the properties of the hyper- 
geometric function (see Appendix) one can easily verify by direct substitution 
into (fll]) that in the range ([5]) the normalized solution of (|T4l) is 

I \ sin (gff) , 1R x 

tti? 1 p {q c -q) 



for < q < q c and 

P(?)=0 (17) 

for q c < q < 1. At that the parameter a and the the fractional power exponent 
(3 must obey the relationship 

a = 2 li (18) 

From here we have 

In a . 
"=te2 (19) 

In Fig. 1 the distribution p(g) is depicted at several values of (3. 



5 Relationship with experimental data 



The main experimental facts on the effect of solvent viscosity on the enzymatic 
reaction rate constant are: 1. the inverse proportionality to fractional power 
exponents of solvent viscosity given by (CQ) and 2. the existence of the limit 
value (3 max of the fractional power exponent with the decrease of cosolvent 
molecular weight given by ([2]). The typical dependence of the fractional power 
exponent (3 on the cosolvent molecular weight M has the form 



^ M ) = W (2°) 
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at M min > M > M max so that (3 max = (3{M min ). The dependence §B) is 
experimentally verified at values of solvent viscosity higher that some minimal 
value {v/Vo)min- For the oxygen escape from hemerythin the experimental 
values are: (3 max = 0.79, b = 1.52, 5 = 0.23, M min = 18, M max = 500000 and 
(r)/m)min ~ 16 ( see M- Substitution of ([20]) into (pl| yields the explicit 
dependence of the parameter a on the cosolvent molecular weight 

a (M, b, 5) = 2 b / MS (21) 

In Fig. 2 the dependence of the parameter a (M, b, 0.23) as the function of the 
logarithm of the cosolvent molecular weight M is plotted for several values 
of the parameter b at 5 = 0.23. In Fig.3 the dependence of the parameter 
a (M, 1.52, 5) as the function of the logarithm of the cosolvent molecular weight 
M is plotted for several values of the parameter 5 at b — 1.52. At M — > M min 
we have a (M) — > a max = a (M min , b, 5). The behavior of a max as the function 
of the parameters b and 5 is plotted in Fig. 4. In the limit a (M, b, 5) — > a max 
the value of the fractional power exponent tends to the limit value (3 max . From 
( fT9~l) we have 

n _ l n a max (00\ 

In 2 

In Fig. 5 the fractional power exponent (3 max given by (f2"2"|) is plotted as the 
function of a max . The experimental value f3 max = 0.79 for the oxygen escape 
from hemerythin from [13] is obtained at a max w 1.75. 



6 Discussion 

A theoretical interpretation of the phenomenon of solvent viscosity dependence 
for the enzymatic reaction rate constant must account for two undisputable 
experimental facts: 1. the inverse proportionality to fractional power exponents 
of solvent viscosity given by ([T]) and 2. the existence of the limit value (3 max of 
the fractional power exponent with the decrease of cosolvent molecular weight 
given by ([2]). Up to now suggested models were mainly concerned with the 
explanation of the first fact while the second one was overlooked or ignored. 
On the other hand the model [IH] suggests the explanation of the second 
fact but encounters difficulties at general description of the first one. The 
present approach is an alternative to all previous attempts that gives combined 
interpretation of both experimental facts in a self-consistent manner. 

The basis of our approach is the relationship (JHJ). The latter is compatible with 
both nonlinear relationship of the viscosity for the movement of the system 
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along the reaction coordinate with solvent viscosity (Q and the linear one 
( 11 II) . The choice among these options is by no means crucial for our approach. 
However we choose the linear relationship (ITTj) because in this case the ap- 
pearance of the fractional dependence (pQ) seems to be more nontrivial than 
in the case of nonlinear relationship. In the latter case we introduce in the 
model the fractional exponent k / 1 from the very beginning by hands. This 
fact obscures the essence of our approach that in fact all additional fractional 
exponents and nonlinearities are superfluous. In (jSJ) we factorize the viscosity 
for the movement of the system along the reaction coordinate C/^o onto the 
contribution from protein matrix and solvent molecules and that from cosol- 
vent molecules. The latter enters in CA/o with the weight q. This weight may 
take different values from the range < q < 1 for the individual samples from 
the ensemble of enzymes. The ensemble is characterized by the distribution 
over the weight p(q). We show that the very existence of this distribution is 
quite sufficient for the appearance of the required dependence of the reaction 
rate constant on solvent viscosity ([I]) consistent with observations. Neither 
the modification of the Kramers' theory nor that of the Stokes law is in fact 
indispensable at such approach. In this regard our approach is conceptually 
much more simple than other interpretations of the phenomenon under con- 
sideration. The approach is motivated by vivid analogy of our problem with 
that of adsorption on heterogeneous surfaces and thus has apparent roots in 
the preceding experience of chemical physics. In this connection we note that 
we can cast (iTTj) into the form 



and call q the " shielding parameter" that shows to what extent protein struc- 
ture shields and screens the reaction coordinate from the solvent viscosity at 
v/Vo > 2. This parameter is less than 1 (0 < q < 1) because in experiment the 
solvent viscosity is varied by the cosolvent of the type of sugars, such as tre- 
halose. It is known that "in trehalose solutions, there is generally a deficiency 
of trehalose and an excess of water in the vicinity of the protein" [46]. The 
deficiency of the cosolvent compared with the bulk in this case is by far takes 
place for protein interior because the larger the cosolvent molecular weight 
so much the worse it can be transmitted by protein structure to the reaction 
coordinate. Then the shielding parameter quantifies the measure of the de- 
ficiency of the cosolvent molecules in the vicinity of the reaction coordinate 
compared with the bulk. We can represent the shielding parameter q in the 
energetic form exp [— ej (/c#T)] where e > 0. In this case instead of the distri- 
bution p(q) over the weight q we obtain the one i/;(e) over the energy e that 
has the form identical to the well known Sips's distribution from the theory 
of adsorption on heterogeneous surfaces. However in our problem the physical 
meaning of the weight q seems to be more lucid than that of the shielding 




(23) 
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parameter. That is why we prefer to deal with the former notion rather than 
with the latter one. 

Our approach takes into account the experimental fact that the fractional 
power exponent [3 tends to some limit value j3 ma x with the decrease of the 
cosolvent molecular weight M. This fact is incorporated in the model as the 
inherent one. The limit value f3 max is the function of the parameters b and S. 
These parameters characterize the protein structure itself and its ability to 
transmit cosolvent molecules to the reaction coordinate. We conclude that the 
limit value of the fractional power exponent f3 max is a unique characteristic of 
the type of the enzyme (that determines the protein structure) rather than 
that of the solution of the enzymes. This conclusion is in agreement with 
experimental data (see Introduction). 

The equation (ITT)) and Fig.l show that the distribution p(q) is zero above 
the critical value of the weight q c = 0.5 given by (IT5l) . We conclude that the 
contribution from cosolvent molecules into the viscosity for the movement of 
the system along the reaction coordinate C/^o can n °t enter with the weight 
greater that one half of that from protein matrix and solvent molecules to 
be compatible with observations expressed by (pQ). For small values of the 
fractional power exponent (3 the distribution is skewed to the upper bound 
q = q c of the range while for large values of j3 it is skewed to the lower bound 
q = 0. In between < q < q c the distribution p(q) tends to infinity at both 
q — > 0+ and q — > q c —. Thus there are two main fractions of the samples in 
the ensemble of enzymes: those with pure contribution from protein matrix 
and solvent molecules (when q « 0) and those with the contribution from 
cosolvent molecules having the weight q w q c and the contribution from protein 
matrix and solvent molecules having the weight « 1 — q c . These fractions 
are the most representative in the ensemble of enzymes. However there are 
always fractions of the samples with intermediate values of the weight q in 
between < q < q c . The behavior of the distribution p(q) in this region is 
quantified by the fractional power exponent (3 (see Fig.l). Thus the physical 
meaning of the fractional power exponent (3 in the experimentally observable 
dependence of the enzymatic reaction rate constant on solvent viscosity ([ID 
is the characteristic of the behavior of the distribution p(q) over the weight 
q in the ensemble of enzymes. This distribution characterizes the solution 
of enzymes, i.e., is determined among with the type of the enzyme by such 
characteristics as the cosolvent molecular weight, pH, temperature, etc. The 
distribution p(q) over the weight q acquires at our approach the status of the 
unique characteristic for the solution of enzymes. 

We conclude that our approach yields conceptually simple interpretation and 
quantitative description of the main experimental data on solvent viscosity 
dependence for enzymatic reactions. 
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7 Appendix 



The formula N2.2.6.15. from [15] is 

a 

J dx x a -\a - xf-\x + z)- p = 



a «+P-± z -PB(a, 13) 2 F 1 (a, p; a + (3; -a/z) (24) 
where | arg z \< 7r; a > 0, Re a > 0, Re [3 > 0. 

The properties of the hypergeometric function used in the paper are: 
see N6.7. in [S] 

2 F 1 (a, b; c; z) = (1 - z)~ a 2 F 1 (a, c - 6; c; (25) 

see N7.2.1.8. in 0Z] 

2 Fi (0, 6; c; *) = 2 F 1 {a, 0; c; z) = 1 (26) 

These properties are used at obtaining (IT6T) . 
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Fig. 1. Distribution p(q) over the weight q (eq. (I16p and eq. (I17p ) at increasing values 
of the fractional power exponent (5: [3 = 0.1 (purple color and the line is skewed to 
the upper bound q = 0.5 of the range); (3 = 0.3 (blue); j3 = 0.5 (violet); j3 = 0.7 
(green); /? = 0.8 (red color and the line is skewed to the lower bound q = of the 
range). 
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Fig. 2. The dependence of the parameter a (M, b, 5) (eq. (|2ip ) as the function of the 
logarithm of the cosolvent molecular weight M at 5 = 0.23 for the values of the 
parameter b from the down line to the upper one: b = 1, b = 1.2, b = 1.4, b = 1.6, 
b= 1.8, 6 = 2. 
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Fig. 3. The dependence of the parameter a (M, b, 5) (eq. (|2T1) ) as the function of the 
logarithm of the cosolvent molecular weight M at b = 1.52 for the values of the 
parameter 5 from the down line to the upper one: 5 = 0.5, S = 0.4, 5 = 0.3, 5 = 0.2, 
5 = 0.1. 
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Fig. 4. The dependence of the parameter draax 

(eq. (|2l|) at M min = 18) as the 

function of the parameters b and 5. 
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Fig. 5. The dependence of the limit value for the fractional power exponent /3, 
(eq. ([22|) ) as the function of the parameter a max . 
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